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Lattice contribution to the electronic self-energy in complex correlated oxides is a fascinating 
subject that has lately stimulated lively discussions. Expectations of electron-phonon self-energy 
effects for simpler materials, such as Pd and Al, have resulted in several misconceptions in strongly 
correlated oxides. Here we analyze a number of arguments claiming that phonons cannot be the 
origin of certain self-energy effects seen in high-T c cuprate superconductors via angle resolved pho- 
toemission experiments (ARPES), including the temperature dependence, doping dependence of the 
renormalization effects, the inter-band scattering in the bilayer systems, and impurity substitution. 
We show that in light of experimental evidences and detailed simulations, these arguments are not 
well founded. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The microscopic pairing mechanism of the high-Tc 
superconductivity remains an unsolved question even 
after twenty years of its discovery. Observations of 
a kink at around 40-70 meV in the band dispersion 
along the diagonal of the Brillouin zone (nodal diec- 
tion) and a peak-dip-hump (PDH) structure at the zone 
boundary by angle-resolved photoemission spectroscopy 
(ARPES^^^^iii^^ have drawn a great deal 
of recent attention as they may shed some light on this 
problem. Although an agreement has been established 
that the kink and PDH structure are signatures of the 
electrons coupled to a sharp bosonic mode, it is still 
strongly debated about the origin of this bosonic mode. 
Influenced by the fact that the high-T c cuprate is a doped 
antifcrromagnctic insulator, a common belief is that this 
bosonic mode has a magnetic 0T\gm££££&iZ&&. How- 
ever, an alternative view is that the electron-phonon cou- 
pling in such a dopcd-insulator can be very strong and 
anomalous because of a number of unusual effects, such 
as poor screening, complex structure as well as the in- 
terplay with correlations. Indeed, oxygen related op- 
tical phonons have been invoked to explain the tem- 
perature and doping dependence of the renormalization 
effect o 10 ' 11 : 12 : 13 : 14 . This idea of phonons being mainly 
responsible for this low energy band renormalization ef- 
fect observed by ARPES has stimulated intense debate. 
There is currently no consensus whether a phonon, a set 
of phonons, or a magnetic mode is the primary cause of 
the band renormalization. 

As mentioned, some important reasons to invoke 
phonon interpretation of the ARPES data are: the 
presence of an universal energy scale in all materials 
at all dopin g 10 ' 15 , particularly in the normal state of 
very low T c materials^; the strong inferred momen- 
tum dependenc e 11 ' 14 ; the existence of multiple bosonic 
mode couplings^ and the decrease in the overall cou- 



pling strength with increased doping, interpreted as a 
screening effect, especially for phonons with eigenvectors 
along the c-axisi 3 -. Yet, there is still a widespread belief 
that phonons are not responsible for the kink features. In 
the following sections, with a comprehensive look at all 
experimental data as well as some recent simulations, we 
address some of the criticisms that have been commonly 
used to argue against the phonon interpretation. These 
include the temperature and doping dependence of the 
renormalization effects, inter-band scattering for bilayer 
system, and the ARPES experiments on impurity substi- 
tuted Bi2212 crystal, Bi 2 Sr 2 Ca(Cu2- a; M 2 ,)0 8+ 5 with M 
= Zn or Ni. Our goal is to clarify these misconceptions 
as being due to oversimplifying the effects of electron- 
phonon coupling in cuprates as well as other strongly 
correlated transition metal oxides. 



II. ASPECTS OF THE ELECTRON-PHONON 
SELF-ENERGY 

A. Temperature Dependence 

In the standard treatment of electron-phonon coupling 
effects, the Debye temperature sets a characteristic tem- 
perature scale, which is well above T c in conventional 
superconducting materials. However in the cuprates and 
other low Fermi energy systems, these two energy scales 
can be comparable. As a result, the temperature depen- 
dence of phonon induced self-energies can be very differ- 
ent from that of conventional superconductors. Accord- 
ing to the ARPES measurements on Bi2212 system, the 
band renormalization in the antinodal region (peak-dip- 
hump structure) shows a dramatic superconductivity- 
induced enhancement when the system goes through a 
phase transition from the normal state to the supercon- 
ducting state. It has been argued that only a mode which 
emerges in the superconducting state and vanishes in 
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the normal state can explain this temperature-dependent 
renormalization effect^i 3 -^. P honoris are thereby ex- 
cluded. 

The sharpness of the renormalization effects due to 
electron-phonon coupling is strongly temperature depen- 
dent, given by the fact that T c of optimally-doped Bi2212 
is close to 100K. To demonstrate this temperature de- 
pendence, we consider the normal state (120 K) and su- 
perconducting state (10 K) of a d-wave superconductor 
coupled to a 36 meV B lg , 55meV oxygen A\ g , and 70 
meV breathing phonon o 14 i 17 . The electron-phonon cou- 
pling for the B\g and breathing phonons are those used 
in Rcf. JA. The A\ g modes involve c-axis motions of 
both planar and apical oxygens, and have two branches 
around 55 and 80 meV. The apical electron-phonon cou- 
pling, derived in Ref. 17 . involves a charge-transfer from 
the apex oxygen into the CuC-2 plane via the Cu 4s or- 
bital, the same pathway that gives rise to bi-layer split- 
ting. However, for simplicity, the apical electron-phonon 
coupling is treated as a constant in the calculations pre- 
sented in this paper. The reason to include three modes 
in the calculation was inspired by the earlier success of 
the two-mode calculation-^ as well as the recent discov- 
ery of multiple mode couplin g 12 ' 13 . For this calculation, 
the tight-binding band structure described in Ref. ,19 
has been used. The real part and imaginary part of the 
self-energy E(k, u>) and the spectral functions A(k,ui) at 
k = (0, 7r) are then obtained within weak coupling Eliash- 
berg formalism^ 4 - and plotted in Fig. [T] Details of the 
calculations are presented in the Appendix. 

At high temperature, both Re£(k, u>) and Im£(k, ui) 
do not exhibit a sharp renormalization feature as shown 
by the dashed curves in Fig. Q](a) and Fig. [lib), respec- 
tively. This demonstrates that the thermal broadening 
effect smears out the sharp renormalization features; in 
addition, broadening effects due to additional many body 
effects would smear out the renormalization features fur- 
ther. Thus, one should not expect to observe any sharp 
renormalization features at k = (0, n) in the normal state 
(~ 100K) from the electron-phonon coupling. In the su- 
perconducting state, the renormalization features of the 
self-energy become much sharper, due to smaller ther- 
mal broadening effect as well as the opening of a su- 
perconducting gap. Consistent with the optimally-doped 
Bi2212 and Bi2223 measuremento 2 ' 4 ' 11 ' 18 , the PDH struc- 
ture of the spectral function at k = (0, 7r) emerges at low 
temperature and disappears at high temperature (nor- 
mal state), as illustrated in Fig. [T^c) and Fig. Hfd), 
respectively. While this behavior is generally expected 
for any phonon, we note that in addition, the self-energy 
from electron-phonon couplings which involve momen- 
tum transfers within and between anti-nodal regions of 
the Fermi surface, such as the apex A\ g and B\ g phonons, 
are greatly enhanced for all k-points due to the large 
density of state enhancements in these regions via the 
opening of a d-wavc gap. A detailed momentum depen- 
dence of the renormalization effects in the normal state 
and superconducting state due to the coupling to the Bi g 
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FIG. 1: The calculated (a) real part, ReE, (b) imaginary 
part of the self-energy, ImE , and the corresponding spectral 
functions, A(k,oj) in (c) normal state and (d) superconducting 
state. An extra 5 meV is added to the imaginary part of the 
self energy in 120K simulation to account for the thermal 
broadening of the quasi-particle life time. The location for 
this calculation is indicated in inset of (a) by a red dot with 
a red curve representing the FS. Insets of (c) and (d) are the 
data of optimally-doped Bi2223 system (T C =110K) taken at 
the normal state (120K) and superconducting state(25K)^, 
respectively. 

phonon has been discussed in Ref. Jd and Ref. \A. Fur- 
thermore, both the dispersion kink and the PDH struc- 
ture in the nodal region have been clearly observed in 
the normal state when measured at a low temperature 
on samples with a lower T^. This lends further support 
to the strongly temperature dependent renormalization 
features due to electron-phonon coupling. 



B. Doping Dependence 

Another problematic statement against the phonon 
scenario stems from the apparent strong doping depen- 
dence of the kink position and strength. Based on the 
wisdom for conventional good metals, phonons should 
not have a strong doping dependence, either in frequency 
of the mode or in strength of the coupling. This is 
not necessary valid for layered, doped insulators with 
strong correlation effects, such as cuprates where dop- 
ing dramatically changes the metallicity and the abil- 
ity of electrons to screen charge fluctuations. We first 
note that many experiments on various cuprates have re- 
ported strongly doping dependent anomalies for several 
phonons, which implies a strongly doping dependent e-p 
coupling for these phonons. For example, from inelastic 
neutron scattering measurements, the breathing mode, 
half-breathing mode, and the bond-stretching modes ex- 
hibit prominent doping dependence of dispersion and en- 
ergy renormalization a 20 ' 21 . In Raman and infrared spec- 
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FIG. 2: The intensity plots of the (a) spectral functions with- 
out resolution convolution and (b) resolution convoluted spec- 
tral functions in the superconducting state (10K) along the 
nodal direction, as indicated in the inset of (b) by the blue 
line. Black curves are the band dispersion extracted from 
the maximum positions of the momentum distribution curves, 
which cut the spectral functions at a fixed energy. The MDC- 
derived dispersions in (a) exhibit three sharp "subkinks" due 
to the coupling to the three phonon modes used in the model; 
while in (b) the subkinks are washed out by the finite instru- 
ment resolution effect leaving an apparent single kink in the 
band dispersion. The white dashed line illustrates the bare 
band for extracting ReE shown in Fig. [3] (a). 



troscopy, the Fano lineshapes of phonon modes in Bi g 
and Bi u symmetry show strong doping dependences^. 
Furthermore, the strength of the phonon energy shift and 
lincwidth variation across T c also changes strongly with 
doping^ 3 -. 

Second, the doping dependence of the renormalization 
effects to the electronic self-energy is sophisticated as in- 
ferred by two recent ARPES studies. One is the observa- 
tion of multiple bosonic mode couplings along the nodal 
direction^. The other is the doping dependence of the 
c-axis screening effect to the coupling between the elec- 
tron and c— axis phonons. As proposed by Meevasana et 
a/i 13 i 24 and Devereaux et al^, for electron-phonon cou- 
pling at long wavelengths, the screening becomes more 
effective at reducing the coupling strength when the c- 
axis conductivity becomes more metallic. Given these 
two results plus the variation of the superconducting gap 
magnitude with doping, the doping dependence of the 
kink energy is highly convoluted in Bi2212 whose super- 
conducting gap has an energy comparable to some of the 
phonons. 

In Fig. [2j we present the intensity plot of calculated 
spectral functions demonstrating a doping dependence 
of the dispersion kink in the superconducting state. The 
superconducting gap size was set to be 40, 20, and 10 
meV for the optimally-doped and more overdoped sys- 
tems. In addition, the coupling strength of the breathing 
mode, whose appreciable coupling is only for short wave- 
lengths and large momentum transfer a 20 ' 21 , remains un- 
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FIG. 3: (a) The ReE extracted from Fig. 2(a) (dashed lines) 
and Fig. 2(b) (solid lines) by subtracting a linear bare band 
(dashed line in Fig. [2ja)) from the band dispersion. The ar- 
rows indicate the maximum positions of the ReE where the 
"single" apparent kink in the band dispersion is usually de- 
fined, (b) Summary of the doping dependence of the apparent 
kink energy and the apparent mode energy extracted by as- 
suming a single mode scenario. 



changed for our doping dependence simulation; while, a 



filter function, uj 2 / (u 



„) with different value of c-axis 



screening frequency uj sc is applied to the c-axis phonons 
(Bi g and Ai g ), to simulate the doping-dependent cou- 
pling strength due to the change of the c-axis screening 
effec t 13 i 24 . We note that although this is a simplifica- 
tion, it represents the general behavior of screening con- 
siderations for phonons involving small in-plane momen- 
tum transfers. Full consideration of screening has been 
given in Ref. and Ref. 24. In addition, a component 

0.005 + uj 2 cV is added in the imaginary part of the self- 
energy to mimic the quasiparticle life time broadening 
due to electron-electron interaction. 

As shown in Fig. [Ha), the coupling to multiple phonon 
modes induces several "subkinks" in the dispersion. The 
positions of these subkinks mostly correspond to the en- 
ergies of phonons plus the maximum d-wavc SC gap, 
An, even through there is no gap along the nodal di- 
rection. This is because when calculating the self-energy, 
one needs to integrate the contributions from the entire 
zone, which makes the electrons in the nodal region " feel" 
the presence of the gap. Furthermore, revealed by the ex- 
tracted real-part of the self-energy, ReE (dashed curves 
in Fig. [3] (a)), the dominant feature in ReE for the OP 
case is induced by 36 meV Bi g mode, while for the OD1 
and OD2 case, the features of the 55 meV Ai g mode and 
70 meV breathing mode gradually out-weight the con- 
tribution from the Bi 9 mode. This demonstrates that 
the change of the SC gap magnitude and the increasing 
screening effect to these phonons because of increased 
doping alters the relative strength of the subkinks caused 
by different modes. 

To simulate the experimental data, we convoluted the 
spectral functions shown in Fig. HJa) with a typical 
ARPES instrumental resolution: 12.5 meV in energy res- 
olution and 0.012 n/a in momentum resolution. As illus- 
trated in Fig. [U^b) and the extracted ReE (solid curves 
in Fig. Ela)), the subkinks are less obvious and become 
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a broadened "single" kink in the dispersion which is lo- 
cated at roughly the energy of the dominant phonon fea- 
ture determined by the maximum position of the ReS 
(the arrows in Fig. EJa)). The doping dependence of 
the kink position is summarized as the solid symbols in 
Fig. [3f b). Assuming a single mode scenario, one can 
obtain the "doping dependence" of the mode energy by 
subtracting out the SC gap size. However, we note that 
this extracted "apparent" mode energy does not match 
any of the modes used in the model; instead, it is a av- 
erage between the dominant features (open symbols in 
Fig. Efb). Clearly, since the kink energy is a sum of the 
superconducting gap and a spectrum of bosonic modes, 
it should not be taken as a precise measurement of the 
energy of any particular bosonic mode. This casts doubts 
to the analysis of the doping dependent properties of the 
kink in the nodal band dispersion based on the single 
bosonic mode coupling scenario^. More importantly, 
this illustrates the complex nature of lattice effects in 
these oxides. 



C. Interband Scattering 

Boriscnko et al££ observed that the scattering rate 
of the bonding and antibonding bands along the nodal 
direction cross each other near the energy of the Van 
Hove singularity, suggesting a strong inter-band scatter- 
ing between the bonding and antibonding bands. They 
argued that only a mode with "odd" symmetry, such as 



spin resonance mode, can mediate such inter-band scat- 
tering. The question whether phonons can induce such 
inter-band scattering has also been raised by these au- 
thors. 

First, we note that recent high energy and momentum 
resolution ARPES experiments on Bi2212 using low en- 
ergy photons( <10 eV) have better resolved the bilaycr 
splitting at the nodal point^. However, as shown in Fig. 
2 of Rcf. ,25, the scattering rate of the bonding and anti- 
bonding band does not exhibit a crossover behavior as 
reported by Borisenko at al.. The inconsistency of the 
data between the two groups implies that more experi- 
ments and better analysis are needed to verify whether 
this inter-band scattering effect is genuine. 

Second, empirically, it has been known for over 15 years 
that interband electron-phonon coupling in the cupratcs 
is very large. The evidence comes from the strong reso- 
nance profiles of many Raman active phonons, which dis- 
play large intensity variations 2 ^. This is generally under- 
stood as a result of strong interband coupling, whereby 
phonons can be brought in and out of resonances via tun- 
ing of the incident photon energy 2 ^. Since, in general, 
phonons can also provide momentum to scatter electrons 
along the c-axis, direct inter bi-layer scattering can occur 
which involves mixing of different symmetries of phonons. 
This can be viewed in a simplified way even if we first 
neglect direct interband scattering and consider a bilayer 
system coupled to c— axis phonons. For q z =0, a simple 
classification of c-axis modes is possible: 



-fc,l,tT L fc,2,cr 



fc,er,a=1.2 



k ■ o 



(1) 



fc,<y,cr,a— 1.2, v 



where a is the index for the electronic states of different 
layers, ei(fc) = 62(h), t±(k) describes the hopping of elec- 
trons between two layers, and the index v can be either 
gerade or ungerade active c-axis modes, with symmetry 
classification with respect to the displacement eigenvec- 
tors to the inversion center of the cell, depicted in Fig. 
SJ After diagonalizing the first two terms by canonical 
transformation, the electron-phonon coupling can be re- 
cast as 



H e-ph = J2 9(9,u)(k,q) { a(g,u)(q)+a\ g u) (-q) 

k,q,(T 



Z k+q,+,cr C k, (+,-), <r C l:+q, -,a°k, (-, + ), a 



h. 



■}(2) 



We have used the c + and c_ for the even and odd 
linear combination of c\ and C2, and subscript g and u 



for the gerade and ungerade mode, respectively. Thus 
for q z — 0, where this classification is possible, ger- 
ade phonons induce intra-band scattering (even chan- 
nel), while the ungerade phonons mediate the inter- 
band scattering process (odd channel) even without di- 
rect electron-phonon coupling across the layers. Yet for 
q z = tt/c, the classification inverts, where gerade modes 
become ungerade and vice- versa, as illustrated in Fig |U 
Thus, even in this simple case, modes at different q z con- 
tribute both to intra and interband scattering, and the 
net weight of the coupling appearing in the self energy is 
then largely determined by the specific momentum struc- 
ture g(k, q). Since the self-energy generally involves sums 
over g 2 , and coupling directly of electrons in adjacent lay- 
ers via phonons are non-negligible, clearly the inter-band 
scattering phenomena can not be used to argue against 
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q z = q z - k/c 
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FIG. 4: The illustration of the gerade and ungerade c— axis 
phonons. The eigenmode of the gerade (ungerade) phonons 
is even (odd) with respect to the mirror plane between two 
CuC>2 layers at q z = 0, while their definition swapped at q z = 
■k/c. The black, grey, and white circles represent the Cu, Ca, 
and O atoms, respectively. 

the phonons being important to the electronic states. 

We also add a remark concerning the electron- 
phonon coupling derived from Raman measurements^ 
in YBa2Cu307 and Bi-2212 compared to that obtained 
from ARPES. While one might naively expect the cou- 
plings to be comparable from Raman and ARPES, we re- 
mark that this situation is remarkably different if the cou- 
pling is strongly moment dependent and whenever corre- 
lations are appreciable. Since Raman measures phonons 
with net zero momentum transfer and ARPES involves a 
sum over all transfers, a sizeable coupling difference may 
be discernablc. This is specifically the case for the B\ g 
phonon, where scattering involving momentum transfers 
across the necks of the Fermi surface near (it, O) 1 ^, further 
enhanced via correlations^, yields a strong contribution 
to the electron self-energy that is absent in phonon self- 
energies. Moreover, a sum rule analysis presented in Ref. 
30 highlights in general how electron and phonon self- 
energies may be qualitatively different in strongly corre- 
lated systems. 

D. ARPES Experiments on Zn and Ni substituted 
Bi2212 

In this section, we comment on recent experiments 
about the renormalization effects in Zn and Ni substi- 
tuted Bi2212 crysta l 31 i 32 . The strength of sharp renor- 
malizations in these substituted crystals is found to be 
weakened compared to the pristine crystals. Since the 
magnetic properties are expectedly modified due to the 
Cu substitution by these impurities, the authors con- 
cluded that the sharp renormalization effects are induced 
by magnetic-related modes, not phonons. 

In fact, a close examination of the data published by V. 
B. Zabolotnyy et al^- and K. Terashima et al^ implies 
that the magnetic property is not the only modification 



due to the substitution by Zn and Ni. First, although 
both sets of data are consistent in the antinodal region 
where the strength of the band renormalization is re- 
duced, they are inconsistent with each other on the kink 
strength along the nodal direction. In the data set of V. 
B. Zabolotnyy et al., the kink strength is weaker in the 
Zn or Ni doped samples, whereas there is no detectable 
change in the data set reported by K. Terashima et al.. 

Second, the data from K. Terashima et al. (Fig. 1(d)- 
(f) in Ref. 32 ) suggest that the bilayer-splitting struc- 
ture is much clearer in the pristine crystals than in the 
Zn and Ni doped crystals. Since the authors have ruled 
out the possibility of a significant doping level difference 
between pristine and impurity-doped crystals, the dis- 
tinct visibility variation of the bilayer structure implies 
a impurity-related change in the electronic structure. 

From these two observations on their data, it implies 
that not only the magnetic properties could change, the 
band structure and scattering behaviors could also be 
affected due to these impurities. It is possible that 
these changes of the electronic structures could "weaken" 
the renormalizaton features observed in the ARPES 
spectrum. Furthermore, we note that the strength of 
electron-phonon coupling could also be modified by the 
substituted impurities: this can be inferred from the 
change of the Fano spectra lineshape of the Bi g 340 
cm -1 phonon in Raman spectral for Zn-doped YBCO^ 3 . 
and Th-doped YBCO 3 ^ resulting from an increase in the 
phonon lincwidth due to impurity scattering. Therefore, 
the experiments on Ni and Zn substituted Bi2212 crys- 
tals are inconclusive experiments to distinguish phonon 
and magnetic modes as the origin of the renormalization 
effects. 



III. CONCLUSION 

We have shown that the temperature and doping de- 
pendence of the renormalization effects, inter-band band 
scattering, and the results of Zn and Ni doped materials 
can be understood in the framework of electron-phonon 
coupling. On the other hand, the issues that make it 
not plausible for the sharp kink being of spin origin, es- 
pecially the spin resonance mode, remain: i) the nearly 
constant energy scale BjS z\ function of doping in small 
gap system 1 ^; ii) the multiple modest; iii) the presence 
of clear kink in the normal state d 13 ' 16 iv) the detailed 
agreement between B lg phonon based explanation of the 
mode coupling as a function of momentu m 1 1 1 14 , while 
the spin resonance with tiny spectral weight (2%) is un- 
likely to give an explanation for both nodal and antin- 
odal renormalization; v) the accumulated evidence for 
lattice polaron effect in underdoped and deeply under- 
doped system o 35 ' 36 . With these weaknesses of the spin 
resonance interpretation, lattice effect is a more plausible 
explanation of the renormalization effects. It remains a 
possibility that the spin-fluctuation and other strong cor- 
relation effects are also very important to determine the 
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electronic structure of cuprates; they likely contribute to 
a smooth renormalization of the band and may be more 
relevant to the higher binding energy. However, opti- 
cal phonons are the most probable origin for the renor- 
malization effects due to sharp modes near 40-70 meV, 
which is also supported by the recent finding of STM 
experiments^^. 
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APPENDIX A: MIGDAL-ELIASHBERG BASED 
APPROACH 

In the calculations presented herein, we evaluate elec- 
tronic self energies and spectral functions via Migdal- 
Eliashbcrg treatment, as discussed in Rcf. 39. The 
dressed Green's function in the superconducting state is 
given in Nambu notation by 



A,, \ = uZ{\j, uj)t + [e(k) + x(K u)]t 3 + 0(k, cj)n 
l ' J [^(k )W )P-[e(k)+x(k,a;)] 2 -^( k)W ) ' I > 

from which the spectral function follows A(k, oj) = — —G'{ 1 (k, ui) as shown in Figs. lc,ld, and 2. The momentum- 
dependent components of the Nambu self energy are given as generalizations of those found in Ref. 39 : 

ujZ 2 (k, uj) = ^Y j M k > P - k)| 2 (K(fy + n f (E p )][5{Lu + Q U - E p ) + S(uj -Q u + E p )} 

+[n b (Sl u ) + n f (-E p )][S(u; - fi„ - E p ) + S(w + tl u + E p )]) (A2) 

X2 (k, u) = '^( k > P - k)| 2 Jr(K(^) +n f (E p )][S(u + fi„ - E p ) - S(lo -(l v + E p )} 

p. V P 

+ [n b (Sl v ) + n f (-E p )][5{cj -fl v - E p ) - 8{lo + ft, + Ep)]) (A3) 

02(k, u) = ^J2 l-9-( k ' P - k)| 2 ^(MO,) + nj{E p )][5(u + Ct v - E p ) - S(lo -tt v + E p )] 

+[n b {n v ) + n f (-E p )][6(u; -Q v - E p ) - 6(u +Q V + E p )]) (A4) 



where v denotes the phonon mode index, and nj and rib 
are the Fermi and Bose occupation factors. <?„(k, q) arc 
the corresponding electron-phonon couplings for mode 
z/, given in reference^ for the B lg and breathing modes. 



I 

We choose to model the A\ g coupling via a momentum 
independent coupling. Further details can be found in 
Ref. 17. 
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